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PREFACE 

This report is one of three reports describing results of investi- 
gations for the "Design, Fabrication/ Exploratory Operation and Opera- 
tional Manual of a Hydraulic Analog Computer". The investigation has 
been sponsored by the Arctic Construction and Frost Effects Laboratory, 
New England Division, Corps of Engineers, U. S. Army, under contract 
No. DA -19-016 -Eng-2743 with the Massachusetts Institute of Technology 
through its Division of Industrial Cooperation under DIC Project No. 
5-7155. 

The entire investigation is directed toward improvement of tech- 
niques for predicting subsurface temperatures, particularly with refer- 
ence to frost and thaw penetration below airfield pavements. This report 
presents the design, construction, and operation of the hydraulic analog 
computer. Two other reports are: (a) Heat Transfer at the Air -Ground 
Interface with Special Reference to Airfield Pavements, and (b) Frost 
Penetration in Multilayer Soil Profiles. 
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INTRODUCTION 


Contract . - 

On 2 October, 1953, the Massachusetts Institute of Technology 
through its Division of Industrial Cooperation entered on a cost-type re- 
search contract with the Arctic Construction and Frost Effects Laboratory, 
New England Division, Corps of Engineers, U. S. Army. This contract. 

No. DA -19-016 -Eng-2743, is entitled "Design, Fabrication, Exploratory 
Operation and Operational Manual of a Hydraulic Analog Computer". 

This is one of three reports covering investigations at M.I.T. during the 
first two years of this contract. 

Specifically, this report describes results of investigations called 
for in the following paragraphs of Appendix A of the contract: 

"1. The project for the first year will consist of making a detailed 
design, preparing working drawings, and constructing a hy- 
draulic analog computer with the following principal features: 

(a) The hydraulic analog computer shall consist of a series 
of small vertical tubes or wells connected to each other at 

the bottom through bundles of capillary tubes or porous stones. 
Each vertical tube shall represent an element of soil at a 
given distance below the ground surface. The surface of a 
liquid in a tube shall represent a temperature if the problem 
is a thermal one, or a head, if a hydraulic problem. The 
capillary tubes or porous stones shall be analogous to the 
thermal or hydraulic resistance to flow. To simulate latent 
heat for freezing and thawing problems, an expansion well 
shall be provided at a temperature level representing the 
freezing point of soil moisture. The size of the well and the 
size of the vertical tubes shall be adjustable to account for 
variable latent heat and volumetric heat. 

(b) The basic hydraulic analog computer shall consist of ap- 
proximately 20 to 30 tubes and shall occupy a space about 3 
feet by 6 feet in area. The computer shall be so designed and 
constructed to permit extension in length without limit, if 
required, and shall be articulated in a series of boards which 
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may be operated either independently or in tandem. 

(c) The hydraulic analog computer shall be a generally appli- 
cable analyzer for one -dimensional heat transfer problems 
both with and without the latent heat condition. 

(d) The hydraulic analog computer shall readily take into ac- 
count variations in soil properties such as changes in thermal 
conductivity, volumetric heat and latent heat with depth in the 
soil, and variation of conductivity and specific heat with soil 
temperature. 

(3) The hydraulic analog computer shall readily take into 
account programmed surface temperature or loading fluctua- 
tions of any sort, including experimentally observed distribu- 
tions. 

2. After completion of construction of the basic computer, the 
contractor shall - 

(a) Make necessary adjustments and modifications to place 
the hydraulic analog computer in satisfactory operating con- 
dition. 

(b) Assemble the hydraulic analog computer at an area at 
Massachusetts Institute of Technology where the equipment 
can be used for the solutions of problems.” 


"5. During the course of the second .years' investigation, the con- 
tractor shall prepare a manual of instructions for future 
operators and equipment, based on the information which has 
been obtained regarding the manipulation of the apparatus in 
the solution of the various problems encountered." 

Scope . - 

The scope of this phase of the contract investigation is, in reality, 
essentially as specified under "Contract". The report is divided into four 
chapters. 

Chapter 1 presents a general discussion of the computer and the 
analogous relationships among parameters in the thermal and fluid flow 
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systems. Further, historical developments in hydraulic analogs both in 
this country and abroad are presented. 

Design requirements and selection of a working fluid comprise 
the major parts of Chapter 2 while the construction and assembly of the 
various components of the computer are discussed in Chapter 3. 

Finally, the programming procedure for the hydraulic analog com- 
puter is presented in Chapter 4. A sample problem solution is given. 


M.I.T. Personnel . - 

Investigations under this contract have been carried out under the 
supervision of Harl P. Aldrich, Jr., Assistant Professor of Soil Mechan- 
ics, Department of Civil and Sanitary Engineering. During the first year, 
Henry M. Paynter, then Assistant Professor of Hydraulic Engineering, 
also supervised parts of the investigation. 

Dr. Ronald F. Scott, Soils Engineer at the Arctic Construction and 
Frost Effects Laboratory, formerly Research Assistant in Civil Engi- 
neering, carried out the majority of the investigation described herein. 

Mr. Rasmus S. Nordal, formerly half-time Research Assistant in Civil 
Engineering, worked on the project during its second year. He contributed 
to the construction and operation of the computer. 

Dr. Aldrich and Dr. Scott, assisted by Mr. George Leung and 
Mr. Nordal, are responsible for the technical writeup. 


Chapter 1 


GENERAL 

The hydraulic analog computer has been designed primarily to 
solve complex one -dimensional thermal problems in soil involving 
freezing and thawing. Of particular interest is the determination of the 
thermal regime below airfield pavements and buildings constructed bn 
non-uniform soil. More specifically, the computer will give the depth of 
freezing or thawing at any time under given soil conditions and boundary 
temperature conditions. 

Before beginning a detailed description of the design, construction 
and operation of the hydraulic analog computer, it will be desirable to 
consider in a general way the manner in which the computer simulates 
the physical conditions of the problems it solves. Further, since other 
analogs have been built in the past, a fairly complete description of their 
construction is also included. The reader will find the remaining parts 
of the report much easier to follow after he reads the following general 
discussion. 

1-01 Description of Computer and Analogous Relationships . - 

Fundamental physical and mathematical principles which are in- 
volved in the one -dimensional thermal problem in soil are well established 
and have been reported elsewhere (1) (2). Mathematical solutions of the 
differential equations have been presented for various simplified cases. 
However, an analytical solution is not practically obtainable for typical 
cases involving non-uniform soil and/or time -varying surface tempera- 
ture. Therefore, these complex problems must be solved by establishing 
a finite difference approximation to the fundamental differential equations 
and solving these numerically or with the aid of an analog. We shall 
describe the latter, in particular the hydraulic analog solution. 

Figure 1 illustrates diagrammatically the analogous relationships 
in the thermal and fluid flow systems while Figure 2 is a photograph of 


"Numbers in parenthesis refer to references in the bibliography. 


the computer. As heat is removed from a given layer of soil its 
temperature falls to the freezing point of soil moisture, normally 32° F. 

The temperature dwells at this point until the latent heat of fusion of the 
soil moisture is removed or, in other words, until the soil moisture 
freezes. After the layer is frozen its temperature drops below freezing. 

A soil profile may be divided into such layers, each layer being repre- 
sented in the computer by a glass tube standpipe with an expansion reser- 
voir. Water level in a standpipe corresponds to temperature for the soil 
layer it represents. Flow of heat within the soil is represented in the 
computing device by the flow of water. The characteristic of each stand- 
pipe and adjustable expansion reservoir, then, is analogous to the thermal 
energy versus temperature characteristic of the soil. 

To represent resistance to heat flow between adjacent soil layers, 
a resistance to fluid flow between standpipes must be provided. In the 
computer, laminar flow takes place in the hydraulic resistance through 
openings of capillary size which can be continuously varied in length. 

Before an actual problem can be solved on the computer, the soil 
profile and its thermal properties must be known. Generally, the thermal 
conductivity, volumetric heat and latent heat of each layer may be esti- 
mated from the water content and dry unit weight of the soil. The computer 
may then be set up with the appropriate standpipes, latent heat reservoir 
openings, and hydraulic resistances. Ground surface temperature, rep- 
resented by the water reservoir to the left in each figure, may be pro- 
grammed automatically to follow actual or assumed variations in surface 
temperature with time by means of a moving paper templet and servo 
motor system. A typical time scale is, for example, the correspondence 
of one minute of computer operation to two days of thermal variations in 
nature. 

The analog computer can handle such refinements as variation of 
volumetric heat and thermal conductivity with temperature. Furthermore, . 
variations in the freezing temperature of soil moisture from layer to layer 
may be easily accounted for by raising or lowering latent heat wells. In 
fact, almost any complex one -dimensional diffusion problem can be solved 
on the computer. Thus consolidation problems in soil engineering and 
diffusion problems in electrical engineering and other fields may be 
studied as well as heat transfer situations. 


The principal advantages of the computer are simplicity of pro- 
gramming and operation, and provision for a continual visual check on 
the solution. While the number of units in the apparatus can be increased 
indefinitely to give more refined solutions, the accuracy of the present 
20 element computer has been clearly demonstrated, as described later, 
by solving problems with known analytic solutions. 

1-02 Historical Developments . - 

Apparently hydraulic analog computers originated almost simul- 
taneously in the United States and in Russia in 1935 and 1936. Moore (3) 
built and patented an apparatus which he called a "Hydrocal" for the 
solution, principally, of heat flow problems through building walls. At 
about the same time Lukyanov (4) built a "Hydraulic Apparatus" in the 
Russian Institute of Roads and Construction which he used to study a 
variety of problems including freezing in soil. A brief description of each 
of the first two analogs is in order, as it will effectively demonstrate the 
principles involved. 

Moore's early Hydrocal was extended and modified in 1939-1940. 
Dr. Scott saw the apparatus during a visit to the University of Michigan 
in 1953. It consisted of eighteen vertical standpipes and connecting re- 
sistances. Moore obtained a variable hydraulic resistance by an appro- 
priate combination, in series and/or in parallel, of calibrated capillary 
tubes. He controlled the temperature of the working fluid which was 
water, by immersing all the resistance tubes in a constant temperature 
water bath. 

Each "standpipe" in Moore's Hydrocal consisted, in reality, of 
two tubes; a fixed vertical manometer tube and a second larger tube which 
could be tilted back to different angles. Thus, the adjustable tube could 
be tilted, within limits, to give an effective cross-sectional area which, 
when combined with the manometer tube, represented the appropriate 
capacitance or volumetric heat. Standpipes were connected to a manifold 
which, with various valves, aided the programming and solution of prob- 
lems. Surface temperature variations were programmed manually by 
raising or lowering a reservoir which moved in guides by pulleys. 

The principal disadvantages of Moore's apparatus were that the 
complex "piping" needed to connect the capillary tubes admitted air 
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which was difficult to remove and rubber connecting tubes perished and 
some corrosion of brass fittings occurred. Moore was aware of the dis- 
advantages of water as a working fluid (meniscus effects, corrosion, 
bacteria growth, etc.) and experimented with other liquids without much 
success. Results of some experiments carried out on the Hydrocal are 
described by Lamb (5). 

The original Russian apparatus described by Lukyanov (4), was a 
good deal broader in its scope as well as in its dimensions. Further, 
Lukyanov states that two other devices were built at the beginning of 
World War II to accommodate two-dimensional flow problems with boundary 
conditions programmed automatically. The writers have found no publica- 
tions which describe these analogs. 

Lukyanov's analog incorporated variable hydraulic resistances 
which apparently consisted of capillary tubes of continuously variable 
length. While the vertical standpipes were of a constant cross-sectional 
area the effective area was varied by inserting prismatic rods. This ap- 
proach leads to the possibility of accounting for variations in volumetric 
heat with temperature by inserting tapered rods. There is no evidence 
that Lukyanov actually did this but the possibility of the refinement is ob- 
vious. 

While Moore's Hydrocal was not designed to account for heat 
stored or released on freezing or thawing, the Russian apparatus did in- 
corporate this latent heat condition. For this purpose an ingenious device 
was attached to each standpipe whereby a hollow cylinder was connected 
to the standpipe by a rubber tube from its base. The cylinder was floated 
in water in such a way that it transferred its liquid volume (latent heat) to 
the standpipe while the fluid level in the standpipe remained constant 
(freezing condition). 

It would appear, from dimensions given in the reference, that 
solutions on Lukyanov's apparatus took a long time. The operating fluid 
was not mentioned specifically but it was presumably water. 

Several hydraulic analogs to diffusion problems have been built 
since 1945. In Sweden, Backstrom (6) (7) studied the cooling of bodies in 
refrigerating cabinets by means of various kinds of simple hydraulic 
analogies using capillary resistance tubes cut to appropriate size. He 
also had occasion to examine the freezing of bodies placed in a refrigerator. 


He used a very simple latent heat device, applicable to one -directional 
heat flow, called a "chicken feed" which carefully meters a given quantity 
of water into the standpipes when the water level corresponding to freez- 
ing is reached. 

In the United States in 1948, Leopold (8) (9) constructed a hydraulic 
analog for studying problems connected with thermal storage, radiation, . 
convection and conduction. This apparatus es located at Pennsylvania State 
University where its use is continued. The resistance tubes are made 
of coiled metal tubing and must be replaced for each problem involving a 
different resistance. This inconvenience is not so serious in Leopold's 
work because the materials through which thermal conduction takes place 
are typical building materials having reasonably well established thermal 
conductivities. Originally the boundary temperature was programmed by 
hand but later a motor driven cam was introduced in order to simulate 
cyclic solar radiation and air temperature. Leopold substituted a sili- 
cone oil (Dow -Corning DC 20) for water as the working fluid. 

In 1948, Barron (10) described a "Viscous Flow Tube Model" which 
was used to study one -dimensional consolidation and swelling in soil. 

He proposed a model design for the freezing or thawing problem where 
each standpipe had an enlargement of suitable dimensions to account for 
latent heat. When water level fell to the enlargement (freezing point) the 
tube was raised a distance equal to the height of the enlargement at a rate 
selected to keep the water elevation constant at the freezing point. 

An air -flow analog computer was built by Coyle (11) in England in 
1951. The air is contained in a closed circuit which also includes a trans- 
former oil, of low viscosity and vapor pressure, which is used to indicate 
the pressure in the air. While air exhibits little variation in viscosity 
with temperature, it has the undesirable properties of having an ex- 
tremely low viscosity and some compressibility. Coyle says the latter 
must be taken into account in accurate work. The resistance tubes are 
metal cupro -nickel and are very small indeed (0.017" bore) because of 
the low air viscosity. They cannot be adjusted continuously. All the 
joints in the analog must be made air-tight and, unlike a water leak, an 
air leak is very difficult to detect. The standpipes were designed to in- 
clude a variation of specific heat with temperature since the problems 
involved heat conduction in metals where the specific heat and thermal 


conductivity may vary appreciably” with temperature. 

At the Massachusetts Institute of Technology, a small hydraulic 
analog was constructed in 1949 (12) (13) to study and demonstrate consoli- 
dation problems in soil. Further, it was used in the settlement analysis 
of the Hayden Library of M.I.T. Again small glass capillary tubes were 
used to simulate the permeability of the soil. 

A paper by Juhasz and Hooper (14) was published in 1953 which 
described a simple hydraulic analog. References to similar analogies 
have been made in articles by de Laclemandiere (15) and Ramachandran 
(16). 


Chapter 2 


DESIGN OF COMPUTER 


This chapter sets forth the basic requirements which must be met 
in the design of the hydraulic analog computer. The selection of a work- 
ing fluid is treated separately because of its important influence on the 
design. The actual design and construction of individual components of 
the computer, together with working drawings, is presented in the fol- 
lowing chapter. 

2-01 Design Requirements . - 

Although the authors gained considerable knowledge of hydraulic 
analog computers from the experience of other investigators, it was 
deemed advisable to make the apparatus as simple and flexible as pos- 
sible because of the research nature of the investigation. This flexibility 
would allow for solutions of other complex diffusion problems in addition 
to freezing and thawing in soil. Further, it was felt that simplicity 
would lead to real usefulness. Consequently, the major design require- 
ments were established as follows: 

(a) The cylindrical standpipe tubes representing volumetric heat, 
should be easily varied in diameter by inserting formers and/ 
or by utilizing tubes of different diameters. 

(b) The hydraulic resistors should be continuously adjustable, 
without having to remove the operating liquid, to accommo- 
date the wide range of thermal conductivity found in soils 
and to enable adjustment during the operation of a problem in 
a case where conductivity varied with temperature. 

(c) An improved and simplified latent heat simulator was neces- 
sary in which freezing and thawing could take place, without 
attention or adjustment of the element, when temperatures 
oscillated about the freezing point. 

(d) If possible, water should not be used as a working fluid be- 
cause of its viscosity variation with temperature and affinity 
for air. A thermostatically controlled temperature bath was 
to be avoided. 


(e) The working fluid should be visible throughout the computer 
so that air bubbles and dirt could be easily detected. The 
components of the computer should, therefore, be transparent. 

(f) Stopcocks for each element must be provided which could be 
operated easily and, if possible, simultaneously (in gangs). 

(g) No part of the apparatus should be susceptible to corrosion 
by the working fluid or outside agencies. 

(h) The fluid meniscus in the standpipe should be eliminated be- 
cause of the difficulty of reading the fluid level and because 
of its small effect on the water movement. 

(i) A programming apparatus which can be operated automatically 
or by hand, should be incorporated for the surface tempera- 
ture variation. 

(j) The computer should consist of about 20 elements or units to 
cover the required range of one -dimensional problems. 

(k) The fotal space required should not exceed 6 or 7 feet by 3 
feet in area. 

Attention was initially concentrated on the resistance tubes which 
are the most important part of the apparatus. However, before the di- 
mensions of the capillary tube through which the operating fluid would 
flow could be established, it was necessary to select a working fluid and 
to determine its physical and chemical properties. 

2-02 Working Fluid . - 

After an extensive investigation of a variety of possible working 
fluids, it was found that water did indeed come closest to fulfilling all the 
desirable requirements. Factors leading to the selection of water are 
outlined in the following paragraphs after the desirable requirements of 
the working fluid have been discussed further. 

Design requirement (b) states that the hydraulic resistor should 
be continuously variable without having to remove and replace the work- 
ing fluid. Thus, once the resistor has been placed in operation the work- 
ing fluid may remain in contact with the resistor for long periods of time, 
until the resistor is removed to make way for another or for cleaning 
purposes. Therefore, the fluid must neither corrode the resistor (or the 
remainder of the apparatus) nor deposit solid or gummy films on the 


resistor passages. Since the dimensions of the passages, even for a 
relatively viscous fluid, must be small to insure 'laminar flow, this places 
severe limitations on both the fluid and material of which the resistor is 
constructed. 

The necessity of excluding air from the apparatus, in particular 
from the hydraulic resistors, placed additional requirements on the 
working fluid. Since the resistor may have a somewhat complicated 
internal shape, it appeared at an early stage that it might be necessary 
to introduce the liquid into the resistor under a vacuum to insure the ab- 
sence of air bubbles. In any event, any simple saturation procedure 
could require large quantities of the fluid, some of which may spill from 
the apparatus. Therefore the fluid should not be expensive or toxic and 
should not leave a residue on evaporation or be gummy or sticky or of- 
fensive in any way. 

Since the apparatus would be exposed to air, although well shielded 
from turbulent contact, the working fluid should not dissolve quantities of 
one of the gases of the air which, under a varying operating pressure 
and temperature, could come out of solution in one of the capillary pas- 
sages and restrict the fluid flow in a manner not allowed for in calcula- 
tions. Dust caps were proposed for the volumetric heat standpipes to 
limit the entry of foreign matter which might also find its way into a re- 
sistor capillary. 

During the solution of a diffusion problem, fluid levels in stand- 
pipes rise and fall and latent heat wells may fill and empty. The working 
fluid must, therefore, drain cleanly from walls of capacitance tubes and 
it must fill and empty the latent heat wells entirely. This implies a non 
surface wetting liquid which is also desirable for reading liquid levels 
accurately because the contact angle between the liquid and wall of the 
standpipes becomes 90°, eliminating a meniscus. 

Since the analog computer will be operated in a laboratory with 
no refined temperature control, it would be desirable to have a liquid 
which has no viscosity variation with temperature over the range of say 
20 to 30°C. Use of such a fluid would eliminate having to immerse the 
hydraulic resistors in a constant temperature bath in order to insure a 
constant resistance during the solution of a diffusion problem. 

Now that the desirable characteristics of a working fluid have 
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been discussed, a list will be given of all the fluids which were considered, 
together with their principal advantages and disadvantages. No extended 
evaluation of each is given, since its possible utility is, in most cases, 
simply disposed of: 

Glycerol: 

Advantages : Does not dissolve air in significant amounts. 

Has a high viscosity and is easily miscible with water to give 
any required viscosity within a large range, which can be 
maintained however, only if the water/ glycerol proportions 
remain constant. 

Disadvantages : Has a high variation of viscosity with tem- 
perature at all water/ glycerol ratios. Glass is wetted by it. 

Carbon Tetrachloride: 

Advantages : Does not entrain air and has a fairly low varia- 
tion of viscosity with temperature. 

Disadvantages: Evaporates readily and is slightly toxic. Has 
a low viscosity and high density. Is a good solvent. 

"Cello solve", "Carbitol" 

and other industrial organic In- 

solvents i* 

Advantages : Have relatively high viscosities which can be ^ 

varied by mixing with water. Do not dissolve and release air. 

Disadvantages : Are strongly solvent and usually have a fairly 
high evaporation rate. They wet glass. 

Kerosene and other fuel oils: 

Advantages : Possess no air -entraining properties and have 
relatively high viscosities. 

Disadvantages : Are inflammable and their viscosity varies 
considerably with temperature. They vaporize easily. 

Silicone oils: 

Advantages : Are stable and non-poisonous. Have a low vari- 
ation of viscosity with temperature. Have a relatively high 
viscosity. 


^Manufactured by Carbide and Carbon Chemicals Company, a Division of 
Union Carbide and Carbon Corporation. 


Disadvantages : Are very expensive, so that any leakage or 
waste would have to be avoided. They wet glass. 

Finally, attention turned back to water which was selected for 
use in the computer. Water dissolves air to a small extent (17) and tends 
to release this air under changed temperature and pressure conditions 
which is not desirable. An investigation was made to see if any additives 
were known which would eliminate this undesirable property. It was 
found that an almost saturated solution of Epsom salts in water did not 
absorb air to the same degree as pure water. However, the difficulties 
incumbent upon using such a solution in the proposed apparatus were 
obvious and were too great to be overcome. Consequently, it was pro- 
posed to use distilled de -aired water and to protect the water from 
turbulent contact with the air. 

Water has a rather low viscosity and thus the dimensions of the 
resistor capillary must be small to insure laminar flow. This in turn 
leads to smaller allowable tolerances in the construction of the resistors. 

Clean water wets clean glass but this disadvantage can be avoided 
by the use of a commercial product called Drifilm.* 

The viscosity of water varies with temperature but to a smaller 
degree than most other fluids, except, the silicone oils. It was 
found that if the computer was placed in a normal laboratory room, away 
from drafts and radiators, it was not likely to undergo large temperature 
changes (of the order of 0.2 to 0.5° C) in the course of an experiment. 

A final precaution was taken to limit the possible growth of bac- 
teria in the water by adding 1/100 of one per cent of mercuric chloride. 
There is no possibility of this quantity of the chemical crystallizing in 
the system and blocking the flow. 

In conclusion, distilled de -aired water has been found to be a 
satisfactory working fluid although it does not possess all of the desirable 
requirements stated originally. 
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Drifilm 9987, General Electric Organosilicon Products. 


Chapter 3 

DESIGN AND CONSTRUCTION OF COMPONENTS 

The hydraulic analog computer, designed and constructed under 
this contract, contains four basic elements, each of which must be de- 
signed and constructed to meet certain requirements. These elements 
are the hydraulic resistors, latent heat wells, volumetric heat stand- 
pipes, and the programming device. Each of these elements will be dis 
cussed in turn after which the connections, stopcocks and frame for the 
computer will be described. 

3-01 Hydraulic Resistors . - 

Two entirely different hydraulic resistors have been used with 
varying degrees of success in the analog computer. The first is shown 
in Figures 5 and 8 and the second in Figure 8. 

a. Design Requirements . - 

The hydraulic resistance R , is defined by the following expres- 

cl 

sion: 

Q - R- ( h i - V 

a 

where: 

Q, rate of flow in cu. in. per min. 

R , hydraulic resistance in min. per sq. in. 

a 

h-^ - h head difference in in. 

In the analog, R must not vary with the head difference or, in other 

a 

words, the rate of flow or velocity must be directly proportional to the 
head difference. This analog requirement, which follows from the pro- 
portionality between heat flow and thermal gradient, implies laminar 
flow. The upper limits of laminar flow in capillaries are discussed in 
Appendix A and experimental procedures for determining the resistance 
R a , are given in Appendix B. 

The principal design requirement involved the selection of a 
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capillary cross-section, which could be easily varied in length, in which 
laminar flow would occur under the anticipated maximum head. Since 
volumetric heat standpipes were to be about 20 inches long, it would first 
appear that resistors should be designed for this condition. In reality 
however, the head difference seldom exceeds 1/2 to 1 inch except for the 
standpipe at the upper boundary in problems involving a step change in 
surface temperature. Even so, in a typical problem having a solution 
time from 5 to 10 hours, a sudden change in head can be adequately rep- 
resented by lowering the boundary reservoir over a period of several 
minutes. Therefore, the resistor capillary was designed to maintain 
laminar flow under head differences of 1 inch. Considerations presented 
in Appendix A indicated that the capillary should have a cross-sectional 
area of between 0.001 and 0.002 sq. in. 

Loss of head at the entrance and exit of the capillary should be 
kept as small as possible. The importance of these losses which are 
proportional to the square of the velocity, depends on their magnitude 
relative to the head loss in laminar flow through the resistor capillary. 

The hydraulic resistors should be constructed if possible, of a 
transparent material so that air could be detected visually. Resistor 
Type A fulfilled this requirement, having been constructed of lucite, but 
resistor Type B was built of stainless steel. 

It was anticipated that each resistor would have a different cali- 
bration curve even though a single calibration for all resistors would be 
desirable. Each resistor should have a graduated scale so that it can be 
easily set to a given resistance. 

b. Hydraulic Resistor Type A . - 

The hydraulic resistor used during the first year of computer 
operation is called Resistor Type A and is shown in Figures 5 and 8. The 
resistance range is approximately 3.0 to 80.0 min. per sq. in. 

Resistor A was designed along the lines of an electrical helical 
resistance. Water flows through a small rectangular space which is 
formed between square threads of a core which is screwed into a barrel. 
In operation, water flows into the resistor chamber at one end, around 
the core in a helical groove until it reaches an exit hole in the side of the 
barrel. The flow is, indeed, reversible. The length of the capillary can 
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be varied by screwing the core within the barrel. The unit is constructed 
of lucite and sealed at the ends by either rubber stoppers or threaded 
lucite caps with 0-ring bushings. 

A graduated scale calibrated in number of turns is attached to the 
barrel. A removable graduated drum made of aluminum serves as a use- 
ful "micrometer" for setting the core position. The graduated drum can 
be readily transferred from one resistor to another. 

The total quantity of water in the resistor remains constant re- 
gardless of the setting. Thus, if the resistance is altered during the 
solution of a problem no sources or sinks are introduced. Before the 
core is rotated the bypass valve must be opened to enable water to pass 
freely from one end of the resistor to the other and thus to prevent a 
pressure difference between the ends. , 

One of the problems in the initial assembly of the resistors is 
their saturation with water. Several methods of accomplishing this were 
attempted. After trying to evacuate all the air from a resistor before 
filling it with water, which proved only moderately successful, one final 
procedure resolved itself into completely immersing both the core and 
the barrel in distilled de -aired water, and assembling the resistor below 
the surface. Another method involved holding the resistor below the level 
of the computer reservoir, which was connected to the bottom of the 
resistor, and allowing water to flow through the resistor continuously 
while assembling it. In this way the presence of water in all the pas- 
sages was assured before the final connections were made. When air did 
occur in a resistor capillary, it was possible to apply an excessive pres- 
sure difference across the resistor, by blowing down a vertical standpipe, 
and to accelerate the water in the capillary to a point where the air bub- 
bles were removed. Any small amount of air which accumulated in parts 
of the resistor where the water was stagnant was not found to be of large 
enough volume to affect the solutions materially. 

The resistance of a pipe of given diameter and length which is 
curved (i.e. has the form of a helix, and through which steady state flow 
is taking place, is different from that of a straight pipe under the same 
conditions. This difference, presumably arising because of the helical 
movement of fluid relative to the pipe itself which leads to a greater 
energy loss, has been investigated by others and described in the 


literature. The importance of this error on the performance of Resistor A 
was investigated. It was found that the ratio of the diameter of the core 
to the diameter of the capillary was too large for a significant difference 
in the friction factor to manifest itself. 

A "pilot" apparatus consisting of six hydraulic resistors of the 
type described was constructed initially. These units were placed in . 
operation and their performance was checked frequently. Two items 
which were of particular interest were the determination of the limiting 
conditions for laminar flow and the general stability of the resistor cali- 
bration. Procedures for determining the hydraulic resistance experimen- 
tally are presented in Appendix B. Since the laminar resistance varies 
directly with the length of the flow path, the calibration curve for the 
resistor should be linear. 

Perhaps the first criterion is to determine the head difference at 
which laminar flow in the resistor breaks down, so that precautions can 
be taken to prevent this situation in a computer solution. If laminar flow 
exists at a reasonably high head under a low resistor setting (short path), 
it follows that flow at higher resistances must also be laminar. There- 
fore, a resistor was tested at two turns. 

For the resistor having a capillary opening of 0.03 in. by 0.05 in., 
it was found that the rate of flow was indeed proportional to the head dif- 
ference in the range in which the apparatus was to operate. A second 
resistor with a capillary opening 0.06 in. by 0.10 in. was built, but on 
testing it was found that it did not exhibit laminar flow at required heads. 

Results of tests showed that the resistance at any one setting did 
not vary with the direction of flow. However, it was discovered that when 
a valve was adjusted to one setting at which its resistance was measured 
and then moved away from it and returned, the new resistance was not 
identical to the previous one. The differences were not large and are 
probably ascribable to the lateral movement of the core of the resistor 
making the capillary path of smaller cross-section on one side and larger 
on the other. 

Further resistance determinations were made after an interval 
of a few days had elapsed, during which the resistor remained at one 
setting, and it was found that no observable change in the resistance had 
occurred. However, later tests have shown small discrepancies in the 
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direction of a slightly larger resistance with time. It was found by con- 
tinual checking of the resistance at one setting, after moving the core 
away from that position and back again, that after a period of time the 
resistance seemed to settle down and remain fairly constant. This may 
be due to a gradual elimination of machining irregularities in the capil- 
lary boundaries or to the gradual formation of an absorbed water film on 
the lucite surface. 

From the series of tests and study of the calibration charts which 
began to accumulate, it was concluded that, for very accurate analog 
solutions it would be necessary to set each resistor experimentally to 
the required resistance, immediately before running the test, and to 
Check it at the end of the test. However, for tests involving rough esti- 
mates of temperature penetration into the soil, it would be sufficient to 
set each resistor according to values obtained from a calibration chart. 

An additional problem encountered with Resistor A was that after 
a long period of time the resistor core began to stick within the barrel 
making the resistance adjustment difficult or impossible. The only prac- 
tical method of remedying this fault, which seemed to be due to small 
distortions or swelling in the lucite, was to take the valve apart, deter- 
mine if possible where the sticking occurred, and remove small shavings 
by machining until the valve operated easily once more. In addition to 
this the valves were greased very lightly with silicone stopcock grease 
before reassembly. In one or two cases the remachining was too exten- 
sive and water apparently moved across the threads in addition to flowing 
in the prescribed path. 

These important disadvantages of Resistor A were the principal 
reasons for seeking a more acceptable design. 

c. Hydraulic Resistor Type B . - 

Hydraulic Resistor B was designed and constructed during the 
second year of computer operation. The final design is shown in Figure 8 
Earlier designs leading to this resistor are described in Reference 20. 
The resistance range for Resistor B is approximately 0.5 to 7 min. per 
sq. in. 

Resistor B consists simply of a hexagonal stainless steel plunger 
moved within a 5/16” O.D. stainless steel tube. Water flows through six 
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capillaries which are varied in length by sliding the plunger within the 
tube. The plunger, which is streamlined at the tip, has a maximum diam- 
eter which is slightly smaller than the inside diameter of the tube. Unlike 
Resistor A, a change in resistor setting does act like a source or sink. 
Thus, a minor error is introduced when the setting is changed during the 
course of a problem solution. 

In order to guide the plunger and prevent rotation within the tube 
which may result in a change of resistance, the hexagonal rod is first 
flattened on 4 of 6 faces so as to fit more precisely the inside of the tube. 
Three pairs of indentations are then made with a center punch into the 
tube wall where the plunger enters it. These provide support and prevent 
rotation of the plunger. 

On calibrating Resistor B it was found that its resistance at 25° C 
was approximately 0.50 min. per sq. in. per inch length of plunger within 
the tube. Figure B-2. The resistor was tested at a low setting (R = 0.54 
min. per sq. in.) and found to exhibit laminar flow to a maximum head h^ 
of 2 inches which is represented in Figure B-2 by a head ratio of about 3. 
These results were considered satisfactory and 10 resistors of this type 
were constructed. 

In the solution of analog problems, it was found that the calibra- 
tion of this resistor did not change appreciably with time. Thus, setting 
the desired resistance from a calibration curve proved to be entirely 
satisfactory. The principal problem appears, in fact, to be related to the 
accuracy with which the length of plunger inserted in the tube can be 
measured. 

At the present time, the analog computer is assembled with the 
first ten resistors of this type. Resistors of Type A are being used, 
where necessary, to give the larger resistances required for thicker soil 
layers at greater depths below ground surface. 

3-02 Latent Heat Wells. - 

A fluid reservoir selected to simulate the stored latent heat of 
fusion of soil moisture in a layer of soil is called the latent heat well. 

The final design adopted for construction is shown in Figures 5 and 7. 
Figure 4 is a photograph showing assembled latent heat wells in operating 
position on the computer. 
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The latent heat well consists of two transparent lucite plates 
separated by a distance of 3/8 inches and centered on the freezing point 
level in the vertical volumetric heat standpipe. The cross-sectional area 
of the well is 0.656 sq. in. The volume of the well is varied by a rectang- 
ular lucite slide which moves between the outer plates. The usable volume 
can be varied from zero to approximately 6 cu. in. 

While this design appears extremely simple, it required consider- 
able experimental testing before an optimum design was selected. Con- 
siderations leading to the design follow. 

It was first decided that none of the latent heat reservoirs de- 
scribed in Chapter 1 would be suitable for the present analog largely 
because of their complexity and lack of flexibility. 

While the soil ideally yields all its latent heat at 32° F when freez- 
ing, in reality freezing may take place over a finite temperature range 
at or below 32° F. Therefore, it was felt that a slight approximation 
could be allowed in the latent heat well in the form of allowing the water 
level to vary slightly while the well emptied or filled. In the final design, 
’'freezing" takes place over a temperature range of about 2°F since a 
typical temperature scale is 0.2 in. equals 1.0°F. The performance of the 
well on draining and filling is shown in Figure 10 while the ideal perform- 
ance and analogous ideal thermal energy diagram for soil are shown in 
Figure 13. 

An important design requirement is that the well must fill com- 
pletely without attention and without trapping air and it must also drain 
freely and completely. In the final design, a tapered groove is provided. 
Figure 7, in the center of the top lucite plate to allow air, which might 
otherwise be trapped at the rear of the well, to escape. To discourage 
water from entering the groove before the air escapes, the slot is coated 
with a silicone non-wetting agent, Drifilm 9987. To aid the exit of water 
from the well when the water level is falling, two grooves were cut along 
the sides of the bottom lucite plate. Short pieces of cotton wick were in- 
serted to hang into the vertical standpipe slightly below the level of the 
latent heat well. These wicks serve to "siphon" water which tends to 
remain along the sides and corners of the reservoir while the water 
drains very slowly. 

At the back of the well, a slide is provided for adjusting the 


capacity of the chamber to accommodate different values of latent heat. 
Since this slide and the well are of a rectangular shape, it was not easy 
to achieve a satisfactory seal at the rear end of the well while permitting 
the slide to be adjusted easily. A solution was obtained by cutting V- 
shaped grooves in the back of the well around the rectangular opening and 
providing matching grooves in a small capping piece of lucite which 
screwed onto the back. A piece of graphite cord was stuffed into the 
groove to serve as a seal. In order to adjust the well setting, it is neces- 
sary to slack the four cap screws slightly, move the slide the required 
amount, and tighten the screws again to maintain the seal. The movable 
slide is provided with a semi-circular notch at the front end. Thus, when 
the slide is closed (zero latent heat) the volumetric heat tube is continuous 
through the freezing temperature. For cases when the area of volumetric 
heat standpipes does not match the area of the semi-circular notch, the 
latent heat well may be removed entirely and a single long standpipe used. 

Latent heat wells are provided for the first 12 elements of the com- 
puter only. The remaining 8 elements are used to establish the lower 
boundary condition in a problem solution. 

3-03 Volumetric Heat Standpipes . - 

Vertical glass tubes about 10 inches long are used to represent 
the volumetric heat capacity of the soil. Tubes of varying diameter 
(typical size, 9/32 in. I.D.), are provided with threaded lucite caps to 
enable the tube to be attached easily to either side of the latent heat wells. 
A seal is accomplished by means of an 0-ring. Loose dust caps, covered 
with wire gauze, are placed on top of the standpipes to prevent contami- 
nation of the water. All standpipes are coated on the inside with Drifilm 
9987, the non-wetting agent, to assist in the reading of the water level 
and to minimize errors due to the meniscus. The volumetric heat stand- 
pipes are shown in Figures 2 and 7. 

In earlier designs, the volumetric heat standpipes were to be 
cemented to the latent heat wells halfway along their 20 inch length. Since 
this required the cutting of a hole and slots in the pipe, it was necessary 
to consider the use of plastic tubing for the purpose. Further, plastic 
tubing was considered to be less breakable. Lucite tubing was thus se- 
lected. 


Lucite tubing of small diameter is extruded and is highly variable 
in diameter which necessitated careful selection to obtain uniform bore. 
When suitable lengths were found they were connected to the latent heat 
wells as indicated above. At this time, and on various other occasions 
during the construction and use of the apparatus, trouble was experienced 
with cracking of the thin lucite tubes, beginning at a stressed point, and 
extending on occasion several inches along the length of the tube in a 
helical fashion. It was discovered that this was partly due to residual 
stresses in the tube due to insufficient annealing before and after ma- 
chining and attachment to screw caps, and partly to solvent migration 
from those sections of the tube which had been cemented to other parts 
of the apparatus. The cracking proved so troublesome that it eventually 
led to the replacement of the plastic tube with glass. 

A manometer tube of small diameter is attached to each volumetric 
heat standpipe for the sole purpose of enabling accurate water level read- 
ings at temperatures around the freezing point. The stored water in this 
manometer tube must be accounted for in the problem scaling as will be 
shown later. 

3-04 The Programmer . - 

The design of a device for programming automatically the reser- 
voir water level to represent any desired surface temperature presented 
one of the most troublesome problems in the computer design. Figure 3 
is a photograph of the programming device and Figures 6, 7, 8 and 9 in- 
clude working drawings of the component parts. 

We seek to transfer a curve of temperature versus time into cor- 
responding water elevations in a reservoir. Since the temperature curve 
may contain steep slopes, a mechanical arrangement involving a stiff 
former with a cam resting on its surface was not considered to be de- 
sirable. Therefore, it was evident from the start that a reversible elec- 
tric motor would be required to operate the reservoir. No difficulty was 
experienced in obtaining a suitable motor and setting up a circuit to oper- 
ate it in either direction. The principal difficulty arose in devising an 
arrangement to inform the motor what to do. 

o After considerable study and several trial designs, it was decided 
that the temperature curve to be programmed would be plotted to an 


appropriate scale on a piece of firm paper such as commercial wrapping 
paper or heavy weight graph paper. The paper would be cut along the 
curve and wrapped around a two drum kymograph which was set at the 
desired speed. The reservoir water level was programmed to follow this 
curve through an electrical probe which was rigidly attached to a float 
indicating the reservoir water level. The probe straddles the paper thus 
interfering with two sets of contact points in the probe. Through two 
small relays, these contacts operate the electric motor to raise the reser- 
voir when no contact is made on either set of terminals, and to lower it 
when both contacts are made. The motor stops when one set of points 
only is in contact. 

It is seen that reservoir water level follows the temperature curve 
by a series of finite steps. The precision of this programming depends 
on the spacing of the contact points, friction in the gearing, inertia of the 
system, and the speed and slope of the paper. When the reservoir level 
was held constant, it was found that as water flows to or from the reser- 
voir it moved about 0.05 in. with each switching operation. 

It is obvious that the reservoir itself does not follow the prescribed 
temperature curve because water level is continually changing relative to 
the reservoir. A reservoir in which water level is maintained at the same 
relative position to it, for example by means of an overflow, would have 
been easier to program. However, it would require a large quantity of 
distilled de -aired water among other things. The computer as constructed 
is a closed system, i.e. no water is removed or added once a problem 
solution begins. 

Because the electrical circuits are such an important part of the 
working of the programmer. Figure 9 has been prepared to show the lay- 
out. Three circuits are involved. A step down from the 115V A.C. supply 
through a transformer and rectifier bridge to a 12.6V R.M.S. full wave 
rectified D.C. supply without smoothing, operates the two small relays 
through the probe contacts. Another circuit involves the relay contacts 
to operate the motor which requires a special circuit of its own to enable 
the reversing to take place. The motor is a 1/20 H.P. split -phase A.C. 
type, with a capacitor start and run, which can be operated in either 
direction by reversing the direction of current through one of the field , 
coils. 


Originally, the rapid switching of the motor, sometimes from 
forward to reverse and back, caused a great deal of trouble by continually 
blowing fuses. This may have been due in part to high transient currents 
which the operations generated. A solution was finally obtained by rede- 
signing the probe, bending the motor relay contacts farther apart to 
minimize the chances of one set making contact before another had 
broken, and inserting small capacitors across the contacts to reduce the 
amount of sparking. It is not known to what extent each of these factors 
contributed to the final elimination of trouble. 

3-05 Connections, Stopcocks, and Frame . - 

Photographs given as Figures 2 and 4 show the component parts 
necessary to connect and control the hydraulic resistors, latent heat 
wells, volumetric heat standpipes, and programmer. Further, the frame 
supporting all components of the computer is indicated. 

Because of the desirability of being able to examine the apparatus 
for air bubbles and since rubber deteriorates or becomes brittle with 
time, "Tygon" transparent flexible plastic tubing of 3/8" inside diameter 
was used for all flexible connections in the apparatus. All the connecting 
tubes were made of this large size to minimize head losses through all 
parts of the system outside the hydraulic resistors. No trouble has been 
experienced with the "Tygon" tubing in computer operations. 

Where it was necessary to control or stop the flow of water be- 
tween two adjacent standpipes and a resistor, a three-way lime glass 
stopcock was used which enabled flow to take place between any two open- 
ings or all three together. To allow the standpipes to be filled at the be- 
ginning of a problem or for calibration purposes, one manifold tube was 
constructed, running the full length of the apparatus. Entry into each 
standpipe is provided by means of an ordinary straight-through lime glass 
stopcock. All the stopcocks have been greased with silicone stopcock 
grease which is a very stable lubricant. 

The components of the computer were first assembled on a frame 
consisting of 1/2 inch diameter aluminum rods connected by screw on 
angle pieces, Figure 4, which gave a very flexible assembly. This proved 
so convenient that it has been used throughout the final apparatus. It 
provides a light-weight, easily movable piece of equipment. New 
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attachments and extra units can be added with very little trouble. 

Since the presence of the latent heat well at the freezing point - 

makes it difficult to read the water level for about 1/2 inch above and be- 
low, an additional tube has been provided beside each standpipe to act as 
a manometer for reading water levels. Behind the manometer is a strip 
of card with two scales: 1) a convenient, frequently used scale containing 
5 Fahrenheit degrees to the inch and having the 32° F mark centered on 
the latent heat well, and 2) an ordinary scale in inches and tenths of an 
inch reading from an arbitrary zero point at the bottom of the standpipe. 

Since a number of these manometer scales were required, it was found 
convenient to have one very carefully and accurately drawn. Photostats 
were made in the number required. The strips have been cemented to 
3/4" wide 0.020" thick aluminum strips. They were then sprayed with 
a clear plastic spray for permanency before mounting. 

The manometer tubes act with the volumetric heat standpipes in 
storing "heat". Therefore, they must be considered as one unit when 
scaling a problem for solution on the computer. 
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Chapter 4 


PROGRAMMING THE COMPUTER 


Programming a problem on the hydraulic analog computer is de- 
fined to include the computational work of scaling the physical problem 
in analog terms, the setup of the computer according to the scaling, and 
the running of the computer to obtain the solution. This chapter will serve 
to describe the programming procedure for the solution of freezing and 
thawing problems in soil. 

4-01 The Scaling Procedure . - 

We seek first to establish relationships between the physical 
conditions of the thermal problem and its analog representation. There 
are a variety of ways in which this scaling procedure may be carried out 
correctly. The most efficient procedure depends on the type of problem, 
available analog equipment, and personal preference of the investigator. 
One procedure will be presented here, in a series of steps to aid the be- 
ginner, for programming a seasonal freezing or thawing problem in 
stratified soil. Reference will be made to a summary of the scaling re- 
lationships, Figure 13, and to the standard scaling sheet. Figure 14. 

These may be folded out for convenient reference. 

Step 1 Data required for the solution are obtained and entered in 
columns 1 to 4 (Figure 14). They are: 

(a) Soil profile with water content and dry density to a depth at 
least equal to the freezing or thawing depth and preferably to 
a depth of 30 to 40 feet or more. The soil profile is plotted 
to scale in column 1. 

(b) Thermal properties of each stratum making up the soil pro- 
file as follows: 

Thermal conductivity, k 
Volumetric heat, C 
Latent heat of fusion, L 

For this example, we shall make no distinction between 
frozen and unfrozen values of k and C. Thermal properties may be ob- 
tained, for example, by procedures given in References 1 and 2. 
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Step 2 The maximum depth of frost or thaw penetration is esti- 
mated. The first 10 to 12 elements of the analog computer are assigned 
to represent this depth since 12 units of the present computer are equipped 
with latent heat wells. In general, the thickness of finite layers of soil, 
called ''lumps", which are represented by elements on the computer, are 
chosen to increase in thickness with depth. Thus, those near the surface, 
where temperature variations and thermal gradients are large, are the 
smallest. 

The thickness Th, of the first soil lump is often chosen together 
with the combined area of the first standpipe and manometer tube. This 
establishes a relationship which must be met between all standpipe areas 
and corresponding lump thicknesses. The volumetric heat for the lump 
is given the notation C where: 

C = C • Th, Btu per °F. (unit area) 


If the lump contains two or more soil types then C is the sum of the prod- 
ucts C • Th’ making up the lump. The thickness of remaining lumps and 
corresponding areas A, must now satisfy the following equation: - 


Th = Th, 
n 1 


c i 

c~ 

n 


A 

n 


or 


— - = Constant 
C 


Standpipe areas are recorded in column 8 for corresponding lumps 
which are numbered consecutively for convenience as shown in column 10,. 
The thickness of lumps is given in column 9 and the distance between 
centers of adjacent lumps is recorded in column 11. 

Step 3 A temperature scale V , is the relationship between the 
prototype temperature and the water elevation, or head, in standpipes: 

V = — , in. per °F. 
s v r 

There are three factors to be taken into account when choosing the tem- 
perature scale. First, the variation in head that simulates the maximum 
temperature variation must be within the height of the standpipes and the 
range of reservoir movement. Second, the scale should be large enough 
to secure adequate precision of the observations. Third, since the 
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selection of V fixes the relationship between thermal energy and volume 
s 

of water, the available volume of the latent heat wells sets an upper limit 

to the temperature scale. A typical temperature scale is 0.2 in. per °F. 

When the temperature scale is chosen, the thermal energy scale 

U , becomes: 
s 

U --zzr ' V , cu. in. per Btu 

s c s 

where the ratio A/C is constant for all lumps. 

The temperature and thermal energy scales are recorded at the 
top of the scaling sheet, Figure 14. 

Step 4 The latent heat of each lump L in Btu (unit area), com- 
puted in the same manner as C, is entered in column 7. The required 
latent heat well volume V^* is then: 

= U g L , cu. in. 

which is entered in column 6. The slide opening required to give this 
volume is obtained from a calibration curve for the latent heat well. 
Figure 11, and recorded in column 5. 

Step 5 A time scale T„, which is the ratio of analog time to 

S 

time in nature, is selected: 

t ■ ' 

T = r— , min. per hr. 

D X 

P 

and recorded at the top of the scaling sheet. A scale which has been 
found convenient is 0.0208 min. per hr. (0.5 min. per day). The time 
scale fixes the value of hydraulic resistance for each resistor, as shown 
in the following steps, and must therefore, be selected accordingly. 

Step 6 The thermal resistance R , between centers of adjacent 

F 

lumps, distance d, is computed from: 

R p = F’ (°F) (hr) per Btu 
which follows from writing Fourier's Law: 
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q = ki = k 


V 1 ‘ V 2 _ V 1 ‘ V 2 
3 R 


If two or more soil types occur between adjacent lump centers, the ther- 
mal resistance is computed for each layer. The sum of these resistances 
equals R^, values of which are recorded in column 12. 

Step 7 The corresponding values of analog resistance R„ , are 
now determined from: 


R a = T s X R p' min ' per sq ‘ in ‘ 

Q 

where the term T g ^ is a constant for the scaling. Values of R & are 
entered in column 13. From calibration curves for resistors, the resistor 
setting is determined for the room temperature at which the solution will 
be run, and is recorded in column 14. 

The theory behind the analog solution, based on a finite difference 
approximation to the differential equations, is more thoroughly treated 
in References 10 and 21. 

Subsurface temperatures vary in accordance with boundary tem- 
perature variations at ground surface and at some selected depth below 
the surface. These boundary conditions must be accounted for in the 
analog representation. 

(a) Upper boundary condition. - The assumed or observed sur- 
face temperature curve constitutes the upper boundary condition. The 
surface temperature programming device, described in an earlier sec- 
tion, is set up to program the end reservoir in the analog to follow the 
required surface temperature curve. 

The surface temperature curve must be drawn to coordinate scales 
in accordance with the temperature scale V g , the time scale T , and the 
available kymograph speeds. The ordinate y in inches becomes: 


y 



• V 


s 


where v g is the surface temperature in degrees F. The abscissa x in 
inches is related to the time in nature t , by: 

x = T • Sp • t 
s P 


< 9 , ■ 




<* 

(« 


where Sp is the kymograph speed in in. per min. Available kymograph 
speeds are given in Figure 12. 

(b) Lower boundary condition . - At a sufficiently great depth be- 
low ground surface the temperature remains nearly constant during the 
year. This depth is known as the level of zero annual amplitude. The 
thermal gradient at this point and below is called the geothermal gradient. 
These factors constitute the lower boundary condition which must be ac- 
counted for in the analog solution. 

The last 8 elements of the analog are provided to account for the 
lower boundary condition. They represent that part of the soil profile 
that is below the zone of freezing. Since the thermal gradients are small 
in this region, we are justified in using larger lumps. Frequently the 
lump spacing is doubled with each increment so it is possible to represent 
adequately a soil profile 50 to 100 feet in depth if necessary. 

The geothermal gradient can readily be taken care of by adjusting 
initially the lower water levels to proper elevations. The capacitances 
of the lower, lumps are so large that little change in water level occurs 
during a solution. If necessary, a metered quantity of water may be added 
to the last reservoir to represent the heat supplied from the earth's 
interior by the geothermal gradient. 

4-02 The Computer Setup and Operation . - 

Once the scaling calculations have been carried out for a problem 
and tabulated, it is necessary to adjust the elements of the computer to, 
the values indicated by the calculations. Although this is not difficult, 
certain techniques have been derived because of the special characteris- 
tics of the computer, to enable the programming to be done as rapidly as 
convenient. The elements to be adjusted are the latent heat wells, the 
volumetric heat standpipes, the hydraulic resistors and the programmer. 

It is generally advisable to set the hydraulic resistors just before the 
problem is run in order to minimize possible changes due to temperature 
errors, etc. 

The correct volume for latent heat wells is set from calibration 
curves. The movable slide is adjusted according to procedures presented 
in Chapter 3. 

Standpipes of appropriate cross-sectional area are selected and 
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incorporated into the computer. If appropriate stan'dpipe sizes are not 
available, it may be possible to insert rods (glass, brass, or stainless 
steel) or other solid forms into the tubes to decrease the cross-sectional 
area to the correct size. 

The surface temperature curve for the programming device is 
drawn on commercial wrapping paper or stiff graph paper, cut out, and 
placed around the kymograph drum. The elevation of the curve in rela- 
tion to the apparatus is adjusted on the kymograph. A final small adjust- 
ment is available in the platform screws which support the entire pro- 
gramming device. The kymograph is set to the desired speed. 

Hydraulic resistors are set on the basis of reliable calibration 
curves, corrected for temperature if necessary. When no curves are 
available, the resistances must be set by trial and error according to 
procedures presented in Appendix B. 

An initial ground temperature profile is assumed and the stand- 
pipes are filled with distilled, de -aired water to heights corresponding 
to the assumed temperatures. The entire apparatus is checked for air 
bubbles and for leaks, both of which must be eliminated. 

Observations made during the analog solution are recorded on an 
Observation Sheet, Figure 15. The entries are obvious. The analog time 
is recorded by means of a stopclock and for tests of long duration a clock 
is used to give the total elapsed time. 

The height of the water level in manometer tubes for volumetric 
heat standpipes is recorded every 5 to 10 minutes depending on the time 
scale and the desired accuracy of the results. In addition, the time when 
each latent heat well starts to drain and the time when' it is completely 
drained, are recorded. These observations are used for determining the 
curve of frost depth versus time. 

Standpipe No. 0 is a piezometer which records the surface tem- 
perature variation which can be compared with that intended to be estab- 
lished by the programmer. 

Finally, the temperature of the water in the apparatus is observed 
at frequent intervals. Precautions must be taken to keep the temperature 
within 0.5 °C of the value for which the resistors are set. 


35 


4-03 Presenting the Computer Results . - 

Observations from the analog computer solution may be plotted in 
a variety of ways, the most common being isotherms on a plot of depth 
versus prototype time. 

The 32° F isotherm only may be plotted to give a curve of depth 
of frost penetration versus time. A direct way of determining the latter 
is based on observations of time when the latent heat well starts to drain 
and when it is completely drained. It is assumed that the lump starts to 
freeze at the former time and is completely frozen at the latter. Figure 16 
illustrates the procedure for determining the frost depth curve. 

Each horizontal line in Figure 16 represents the beginning and end 
of freezing for the respective lump. If a latent heat well is not completely 
drained before the adjacent one begins, half of the overlap time is in- 
cluded in the time for drainage of both wells. If different amounts of 
latent heat are represented by the two wells, the overlap time is divided 
in proportion to the volume of the wells based on the assumption that 
freezing in this region is proportional to time. This overlap, which is 
theoretically not possible, occurs because the latent heat will discharge 
its volume under a slightly varying head. 

When a well is partly drained, the depth of freezing in the soil 

a 

lump is assumed to be proportional to the amount of water that is drained 
from the well. 

4-04 Verification of Computer Results . - 

Operation of analog computers of any type are subject to errors 
which may be contributed solely by the operators, such as errors in pro- 
gramming or which may be caused by mechanical difficulties in the com- 
puter, such as leakage and blocking of resistors by air bubbles or 
impurities. Also, an analog solution obtained by ’'lumping" time or space 
variables involves some degree of inaccuracy, which is inherent with the 
finite difference approximation. 

One may achieve a verification of computer results by reprogram- 
ming the problem to different scales. If the results turn out to be in 
agreement with the first set, there is good indication that the results are 
correct. This procedure is particularly favorable in checking errors in 


scaling and setting. However, it involves time and labor. Probably satis- 
factory results can be obtained by simply double checking all initial scal- 
ings and settings before and after the problem is run. 

The foremost method for checking errors caused by mechanical 
difficulties is by constant calibration of all contributing components of the 
computer. Any mechanical difficulty will be noted immediately if deviation 
is observed in the calibration chart. 

Indeed, analog results can always be verified by comparison with 
results of analytical or numerical solutions. The drawback is that ana- 
lytical solutions exist for only a limited number of heat flow problems, 
and thus for complicated soil profiles and boundary conditions, numerical 
solutions will involve tedious labor. Indeed, it is because of such draw- 
back that the hydraulic analog computer is designed. 

The analog computer inherits an inaccuracy in representing soil 
in lumps. The degree of such inaccuracy depends among other things on 
the size of each lump and the magnitude of thermal gradients. A study 
of the magnitude of this inaccuracy was carried out by Scott (21) by com- 
paring analog results of different lump sizes with results of an analytic 
solution, and it was found that the error involved is negligible when the 
solutions are programmed according to the procedure recommended in 
Chapter 4. 


In summary, the hydraulic analog computer, the design and opera- 
tion of which is presented in this report, has been found to be a very useful 
tool for solving complex one -dimensional freezing and thawing problems 
in soil. Indeed, its applicability for solving other diffusion problems in 
engineering is practically limitless. Thus, to use the computer to solve 
consolidation problems in soil, for example, one needs knowledge only of 
the analogous physical relationships. The programming steps are similar. 
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PROBe CONTACT I 



a. PROBE AND RELAY CIRCUITS 



b. MOTOR CIRCUIT 


Symbol 



D 

D 1 ' D 2 /D 3 
E 


F 1 ,F 2 
G 

°i ,g 2 

L 1 ' L 2 
M 1 , M 2 

P 


Electrical Elements Key 

! 

1 Element 

i 

Field coils of 1/20 h.p. A.G . split-phase capacitor start and run motor. 

i 

. j 

5 mf. motor capacitor. 

3PDT motor relay coil. ! 

i 

I 

Contacts of motor relay. 

0.002 mf. capacitors across motor relay contacts to inhibit sparking. 
First probe relay coil. 

Contacts of first probe relay! 

j . 

Second probe relay coil. i 

Contacts of second probe relay. 

Limit switches to prevent reservoir over-running. 

i 

Up and down switches for motor, operated by probe relays. 

115V. A.C. power supply. 


A 5A. fuse is also inserted in the! motor circuit. 

j ■ . ' 

Two push-button switches are used to control the starting position of the 
reservoir. 1 


h 

I 
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VOLUME IN CU. IN. 


FIGURE 10 


VOLUME IN CU. INCHES 




PERIMETER SPEED IN INCHES PER SECOND 



a. IDEAL STORED WATER VS HEAD I 
DIAGRAM FOR AN ELEMENT OF j 
THE COMPUTER REPRESENTING ! 
SOIL "LUMP". 


Head, h 
in inches 


above latent- 
heat well 


below latent 
heat well 



Stored Water, V 
in cu . in ■ 


A, Standpipe and manometer area, sq .in 

\ 

Vjy Latent heat well vol cu.in. ' 
R q , Analog resistance, min .per sq . in . 


t , Analog time, min. 


b. IDEAL THERMAL ENERGY VS 
TEMPERATURE DIAGRAM FOR A 
SOIL "LUMP" OF UNIT AREA 
(SQ.FT.) AND THICKNESS Th. 


Temperature, v k 
in 0 F * 



Thermal Energy, u 
in Btu 


C, Volumetric heat, Btu per°F 

L, Latent heat, Btu 

R q , Thermal resistance, °F hr. per Btu 

t , Time in nature, hr. 
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SCALING RELATIONSHIPS 
ANALOG -r TO — PROTOTYPE 


I b 

Temperature Scale, V ; in. per °F 


Thermal Energy Scale, U s = — • V ; cu . in. per Btu 


t : R A 

Time Scale,' T = — =) — ; min . per hr. 
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FIGURE 13 




HYDRAULIC ANALOGUE COMPUTER 

SCALING SHEET TceT q 

I TEST NO. : O 

PROBLEM = f rro ^ Penetration be/or/ /)sf>ha/t Concrete Pa re n erf date : 6 Jo/y, /45S 
f° r Ctcp C horge /n Surface ff/rperafare computed by: (Z/o 


TIMESCALE I MIN = __^_ DAYS 
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PROFILE 
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Haai 
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28 

850 

23 
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s p ™ e d ‘ ANALOGUE R p r * Stih’ng 
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BTU | (INCH*) PROFILE 6 TU IN 2 ( / HO i £ 2 D 
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2/3 .oS& 


2(3 .056 


0.44 


/.33 !.o45 53~J .10*1 


/Zoo Z5 /.30 A 33 
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43/ 2 83 (45o .z i5 


4.3 1 


4-31 


5.(4 


ZZS 


Q2ZS 


t]3o . 2^7 











8. 1 


/ 3.8 




\mm\ 


2 \ /3 \ /<*■ 


.ISO I 0.34 I /& '8 



•s , 

.So 


N .15 0.64 , 5.J8 


0.64 >5.30 


.544 0.14 ,5?.C, 


.344 o.~J4 i5.4l 


■344 0 .j4 iS.qc, 


0-V 1551 


/oq )5.4o 


1.00 iSSZ 


(Taels') 
( ( S C.3S 


f.ze c.3-j 


/ 63 O. 6o 


(56 o.5o 


363 i.5o 


56o 455 


6.45 Z50 



n . 60 


5 -7 6 


^ -IlS 


^ 1-7° 


$ 2-63 


\ J-35 


^ (3 fa 


iS-JO Q.oo 
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5Z.Z ZZ .00 



FIGURE 14 


tfes/sToa. # — : i?£5/s to £. 









































































































































FIGURE 15 
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HYDRAULIC ANALOGUE COMPUTER observed by: test ho.JL 


Frost Pe"£re4 r /°N sbloia/ 0sph*cj ^v < r2^# ERVAT,0N SHEET oate sheet_L_of_L_ 

PROBLEM • f^^ateMT FQ£ s r£P chas/££ f* SoeFAce Jt/ypeeAT-oee 


DATE 

AND 

HOUR 

ELAPSED TIME 

J— MIN> ^_0*Y8 

STANDPIPES : TEMPERATURE SCALE INCH = 1° F 

REMARKS 
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m 
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Appendix A 


UPPER LIMITS OF LAMINAR FLOW'IN CAPILLARIES 

From considerations of laminar flow through openings of any 
geometric cross-section (18): 

V R „ 2 

Q = c w g - i a (A-l) 

S |JL 

where: 

Q, rate of flow in cu. cm. per sec. 
y , unit weight of water in g. per cu. cm. 

Rr* hydraulic radius in cm. 

|jl, viscosity of water in g. sec. per sq. cm. 
i, hydraulic gradient (dimensionless) 
a, cross-section area in sq. cm. 

and 

C , a shape constant 
s 

= 1/2 for a circular capillary tube 
= 1/3 for closely spaced parallel plates. 

In tubes of circular cross-section, laminar flow will occur if: 

v d v 

* W 

Reynolds number, R = < 2000 - - - - (A -2) 

® M* § 

where: 

v, velocity in cm. per sec. 
d, tube diameter in cm. 

At the upper limit of laminar flow, R e = 2000, the relationship 
between hydraulic gradient and tube diameter becomes: 


A-l 


id = ™ii l — t 
s 


(v*h) 


- - - - (A-3) 


Since: 


then: 


Rtt = x and C = lr 
H 4 s 2 


, ,3 


id = 4000 g ( -t- 


m 


-5 

At a temperature of 70°F, jx = 10 and therefore: 

id 3 = 6.31 x 10~ 3 cm. 3 


= 3.84x 10 ' 4 in. 3 


(A -4) 


For closely spaced parallel plates the upper limit of laminar flow 
(19) occurs when R g = 1000 where d is defined as the distance between 
plates. In this case: 

id 3 = 0.722 x 10' 4 in. 3 - - - - (A-5) 

Curves of i vs d are given in the accompanying Figure A-l. A 
curve for a rectangular cross-section will lie somewhere between the two 
curves, approaching that for circular tubes when the cross-section is 
square. These curves will serve as rough guides for determining limit- 
ing dimensions for resistor capillaries. 
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DIMENSION d IN INCHES 


RELATIONSHIP BETWEEN 
HYDRAULIC GRADIENT AND DIMENSION d 

FOR UPPER LIMIT OF LAMINAR FLOW THROUGH 
CIRCULAR TUBES AND BETWEEN PARALLEL 
PLATES FOR WATER AT 70°F 

FIGURE A - I 


Appendix B 


CALIBRATION OF HYDRAULIC RESISTORS 


The calibration of a hydraulic resistor may involve an experi- 
mental determination of the resistor setting required to give a desired 
value of resistance R . If a calibration curve is desired, the experimen- 

cl 

tal procedure seeks to determine the resistance at several resistor 
settings so that a curve may be drawn, such as that shown in Figure B-2a. 

A number of experimental procedures are available for calibrating 
a resistor. The most useful procedure depends on the investigator, the 
objective of the calibration and the physical circumstances under which 
the resistor exists. 

In all cases the temperature of the water must be measured at the 
time of the calibration and all resistance determinations which are to be 
compared must be reduced to a common temperature by the following 
equation: 

R al = hL_ • 

^a2 ^2 

where p is the viscosity of water. 

B-01 Constant Head Test. - 

A test setup for determining the hydraulic resistance by the con- 
stant head method is shown diagrammatically in Figure B-la. The pro-*- 
cedure consists simply of setting the resistor at the desired position, 
raising the headwater reservoir to give the required head loss h^ - h^, 
and determining the rate Q, at which water flows through the resistor. 
Thus: 



It is generally felt that one of the following procedures, which 
make use of a falling head, represents a more desirable approach. 


B-l 


B-2 


B-02 Falling Head Tests. - 


Two falling head test procedures for calibrating hydraulic resistors 
are described below. These have the advantage that the resistor can be 
calibrated "in place" in the computer making use of a volumetric heat 
standpipe already connected to the resistor. In both cases the general 
setup for the test is shown in Figure B-lb. The resistor is connected to 
a standpipe on one side and to a reservoir which is suspended on a cali- 
brated spring on the other. The reservoir has a cross-sectional area such 
that a rise in the water level produces an elongation of the spring of the 
same magnitude so that the water level remains at the same absolute 
elevation. A graduated scale is attached to the standpipe to facilitate the 
falling head reading. As an alternate to using the reservoir and spring 
combination, the tailwater reservoir can be any container such as the 
programming reservoir in the computer which can be easily adjusted to 
keep the tailwater elevation constant. 

Water flows through the resistor, when the stopcock is opened, 
under a variable head h. We may write the following differential equation 
by continuity: 




- - - - (B-l) 


where: 

A, cross-sectional area of standpipe 
t, time 


Rearranging and integrating for the limits: 


and solving for R a : 


t = 0 when h = h 

o 

t = t^ when h = h^ 



A In 




(B-2) 


R = 
a 


'1 


2.3 log ° 

*1 


(B-3) 


or: 


B-3 


4 


ft 


(a) Falling head test with single observation . - It is frequently 
convenient to let 

h 

h l ' 2 

in which case Equation (B-3) becomes: 


R = 1.445 -£ 
a A 


While this procedure is very quick, the result depends on one 
observation only. Indeed, it is possible to repeat the test but the follow- 
ing procedure offers additional advantages. 

(b) Falling head test with multiple observations . - Equation (B-3) 
shows that if R & is constant, t^ is proportional to log h o /hj . It follows 
then that a plot of t^ to a natural scale vs. h / h^ to a logarithmic scale 
will yield a straight line for laminar flow. 

This latter fact is utilized where the falling head is recorded as a 
function of elapsed time. The observations are plotted on semi -logarithmic 
paper and the resistance is computed from Equation (B-3). The term 



is the time intercept over one log cycle. An example is shown in Fig- 
ure (B-2b). 

This procedure is more elaborate than that based on a single ob- 
servation of time. However, since the result is based on a series of 
observations it is far more dependable than the former. Moreover, the 
test can be used to check for laminar flow through the resistor. For a 
given capillary length there exists a critical gradient above which the 
laminar nature of the flow breaks down and becomes turbulent. In this 
case the semi -log curve deviates from a straight line and R is no longer 
independent of head. 

Since the resistors normally operate at heads of less than one inch, 
the calibration should be carried out at heads of no more than 3 or 4 inches 
in order to minimize any small dependence of resistance on head. 
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